function mpc = case9
%CASE9    Power flow data for 9 bus, 3 generator case.
%   Please see CASEFORMAT for details on the case file format.
%
%   Based on data from Joe H. Chows book, p. 70.

%   MATPOWER

%% MATPOWER Case Format : Version 2
mpc.version = '2';

%%-----  Power Flow Data  -----%%
%% system MVA base
mpc.baseMVA = 100;

%% bus data
%	bus_i	type	Pd	Qd	Gs	Bs	area	Vm	Va	baseKV	zone	Vmax	Vmin
mpc.bus = [
	1	3	0	0	0	0	1	1	0	345	1	1.1	0.9;
	2	2	0	0	0	0	1	1	0	345	1	1.1	0.9;
	3	2	0	0	0	0	1	1	0	345	1	1.1	0.9;
	4	1	0	0	0	0	1	1	0	345	1	1.1	0.9;
	5	1	90	30	0	0	1	1	0	345	1	1.1	0.9;
	6	1	0	0	0	0	1	1	0	345	1	1.1	0.9;
	7	1	100	35	0	0	1	1	0	345	1	1.1	0.9;
	8	1	0	0	0	0	1	1	0	345	1	1.1	0.9;
	9	1	125	50	0	0	1	1	0	345	1	1.1	0.9;
];

%% generator data
%	bus	Pg	Qg	Qmax	Qmin	Vg	mBase	status	Pmax	Pmin	Pc1	Pc2	Qc1min	Qc1max	Qc2min	Qc2max	ramp_agc	ramp_10	ramp_30	ramp_q	apf
mpc.gen = [
	1	0		0	300	-300	1	100	1	250	10	0	0	0	0	0	0	0	0	0	0	0;
	2	163	0	300	-300	1	100	1	300	10	0	0	0	0	0	0	0	0	0	0	0;
	3	85	0	300	-300	1	100	1	270	10	0	0	0	0	0	0	0	0	0	0	0;
];

%% branch data
%	fbus	tbus	r	x	b	rateA	rateB	rateC	ratio	angle	status	angmin	angmax
mpc.branch = [
	1	4	0	0.0576	0	250	250	250	0	0	1	-360	360;
	4	5	0.017	0.092	0.158	250	250	250	0	0	1	-360	360;
	5	6	0.039	0.17	0.358	150	150	150	0	0	1	-360	360;
	3	6	0	0.0586	0	300	300	300	0	0	1	-360	360;
	6	7	0.0119	0.1008	0.209	150	150	150	0	0	1	-360	360;
	7	8	0.0085	0.072	0.149	250	250	250	0	0	1	-360	360;
	8	2	0	0.0625	0	250	250	250	0	0	1	-360	360;
	8	9	0.032	0.161	0.306	250	250	250	0	0	1	-360	360;
	9	4	0.01	0.085	0.176	250	250	250	0	0	1	-360	360;
];

%%-----  OPF Data  -----%%
%% generator cost data
%	1	startup	shutdown	n	x1	y1	...	xn	yn
%	2	startup	shutdown	n	c(n-1)	...	c0
mpc.gencost = [
2	1500	0	3	0	0.0 	0; % modified
2	1500	0	3	0	0.0 	0; % modified
2	1500	0	3	0	0.0 	0; % modified
];


%% candidate dc bus data
%column_names%   busdc_i grid    Pdc     Vdc     basekVdc    Vdcmax  Vdcmin  Cdc
mpc.busdc_ne = [
		1              1       0       1       345         1.1     0.9     0;
		2              1       0       1       345         1.1     0.9     0;
    3              1       0       1       345         1.1     0.9     0;
    4              1       0       1       345         1.1     0.9     0;
	  5              1       0       1       345         1.1     0.9     0;
		6              1       0       1       345         1.1     0.9     0;
		7              1       0       1       345         1.1     0.9     0;
		8              1       0       1       345         1.1     0.9     0;
		9              1       0       1       345         1.1     0.9     0;
];


%% candidate branches
%column_names%   fbusdc  tbusdc  r      l        c   rateA   rateB   rateC status cost
mpc.branchdc_ne = [
1	4	0.0	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;
4	5	0.017	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;
5	6	0.039	0.00	 0.00  150	 0.0	 0.0	 1.0	 1;
3	6	0.0	0.00	 0.00  300	 0.0	 0.0	 1.0	 1;
6	7	0.0119	0.00	 0.00  150	 0.0	 0.0	 1.0	 1;
7	8	0.0085	0.00 	 0.00  250	 0.0	 0.0	 1.0	 1;
8	2	0.0	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;
8	9	0.032	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;
9	4	0.0	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;

1	4	0.0	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;
4	5	0.017	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;
5	6	0.039	0.00	 0.00  150	 0.0	 0.0	 1.0	 1;
3	6	0.0	0.00	 0.00  300	 0.0	 0.0	 1.0	 1;
6	7	0.0119	0.00	 0.00  150	 0.0	 0.0	 1.0	 1;
7	8	0.0085	0.00 	 0.00  250	 0.0	 0.0	 1.0	 1;
8	2	0.0	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;
8	9	0.032	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;
9	4	0.0	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;

1	4	0.0	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;
4	5	0.017	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;
5	6	0.039	0.00	 0.00  150	 0.0	 0.0	 1.0	 1;
3	6	0.0	0.00	 0.00  300	 0.0	 0.0	 1.0	 1;
6	7	0.0119	0.00	 0.00  150	 0.0	 0.0	 1.0	 1;
7	8	0.0085	0.00 	 0.00  250	 0.0	 0.0	 1.0	 1;
8	2	0.0	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;
8	9	0.032	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;
9	4	0.0	0.00	 0.00  250	 0.0	 0.0	 1.0	 1;
];

%% candidate converters
%column_names%   busdc_i busac_i type_dc type_ac P_g   Q_g  islcc  Vtar rtf xtf  transformer tm   bf filter    rc      xc  reactor   basekVac Vmmax   Vmmin   Imax    status   LossA LossB  LossCrec LossCinv  droop Pdcset    Vdcset  dVdcset Pacmax Pacmin Qacmax Qacmin cost
	mpc.convdc_ne = [
		1       1   1       1       -360    -1.66   		0 1.0        0.001  0.1 1 1 0.08 1 0.001   0.09 1  345         1.1     0.9    100     1     	1.1033 0.887  2.885    2.885       0.0050    -52.7   1.0079   0  1000 -1000 500 -500 3;
		2       2   1       1       -360    -1.66   		0 1.0        0.001  0.1 1 1 0.08 1 0.001   0.09 1  345         1.1     0.9    100     1     	1.1033 0.887  2.885    2.885       0.0050    -52.7   1.0079   0  1000 -1000 500 -500 3;
		3       3   1       1       -360    -1.66   		0 1.0        0.001  0.1 1 1 0.08 1 0.001   0.09 1  345         1.1     0.9    100     1     	1.1033 0.887  2.885    2.885       0.0050    -52.7   1.0079   0  1000 -1000 500 -500 3;
		4       4   1       1       -360    -1.66   		0 1.0        0.001  0.1 1 1 0.08 1 0.001   0.09 1  345         1.1     0.9    100     1     	1.1033 0.887  2.885    2.885       0.0050    -52.7   1.0079   0  1000 -1000 500 -500 3;
		5       5   1       1       -360    -1.66   		0 1.0        0.001  0.1 1 1 0.08 1 0.001   0.09 1  345         1.1     0.9    100     1     	1.1033 0.887  2.885    2.885       0.0050    -52.7   1.0079   0  1000 -1000 500 -500 3;
		6       6   2       1       -360    -1.66   		0 1.0        0.001  0.1 1 1 0.08 1 0.001   0.09 1  345         1.1     0.9    100     1     	1.1033 0.887  2.885    2.885       0.0050    -52.7   1.0079   0  1000 -1000 500 -500 3;
		7       7   1       1       -360    -1.66   		0 1.0        0.001  0.1 1 1 0.08 1 0.001   0.09 1  345         1.1     0.9    100     1     	1.1033 0.887  2.885    2.885       0.0050    -52.7   1.0079   0  1000 -1000 500 -500 3;
		8       8   1       1       -360    -1.66   		0 1.0        0.001  0.1 1 1 0.08 1 0.001   0.09 1  345         1.1     0.9    100     1     	1.1033 0.887  2.885    2.885       0.0050    -52.7   1.0079   0  1000 -1000 500 -500 3;
		9       9   2       1       -360    -1.66   		0 1.0        0.001  0.1 1 1 0.08 1 0.001   0.09 1  345         1.1     0.9    100     1     	1.1033 0.887  2.885    2.885       0.0050    -52.7   1.0079   0  1000 -1000 500 -500 3;
];
